Transition metal dimers as potential molecular magnets: A challenge to 

computational chemistry 
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Dimers are the smallest chemical objects that show magnetic anisotropy. We focus on 3d and 
4d transition metal dimers that have magnetic ground states in most cases. Some of these mag- 
netic dimers have a considerable barrier against re-orientation of their magnetization, the so-called 
magnetic anisotropy energy, MAE. The height of this barrier is important for technological appli- 
cations, as it determines, e.g., the stability of information stored in magnetic memory devices. It 
can be estimated by means of relativistic density functional calculations. Our approach is based on 
a full-potential local-orbital method (FPLO) in a four-component Dirac-Kohn-Sham implementa- 
tion. Orbital polarization corrections to the local density approximation are employed. They are 
discussed in the broader context of orbital dependent density functionals. Ground state properties 
(spin multiplicity, bond length, harmonic vibrational frequency, spin- and orbital magnetic moment, 
and MAE) of the 3d and 4d transition metal dimers are evaluated and compared with available ex- 
perimental and theoretical data. We find exceptionally high values of MAE, close to 0.2 eV, for four 
particular dimers: Fe2, C02, Ni2, and Rh2. 

PACS numbers: 31.15.E-, 31.15.aj, 75.30.Gw, 75.50.Xx 

Keywords: density functional calculations, 3d and Ad transition metal dimers, magnetic anisotropy energy, 
MAE, orbital polarization correction, OPC 
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I. INTRODUCTION 

During the past decades we have been eyewitnesses of 
enormous progress in experimental material design. It 
has become common that this experimental progress is 
accompanied by material-specific theoretical description. 
Beyond this, theoretical modelling has even partly re- 
placed expansive experimental scannings of large series 
of samples with, e.g., different compositions. This ad- 
vance has become possible due to the large predictive 
power of quantum mechanical calculations, in particular 
of density functional approaches. 

Of particular scientific and practical interest is the 
study of small particles. In the realm of chemistry, the 
smallest considered objects are dimers. While many main 
group elements form natural homonuclear dimers at am- 
bient conditions, the production and experimental study 
of metal dimers is difficult. On the other hand, as will 
be demonstrated in the present work, some metal dimers 
show exceptional magnetic properties that are of poten- 
tial interest for, e.g., magnetic information storage tech- 
nology. It is a major challenge to computational chem- 
istry to identify the most promising homo- and heteronu- 
clear dimers that deserve intense experimental study and 
further theoretical modelling, e.g., related to a possible 
embedding into an appropriate matrix. In the present 
study, we will concentrate on homonuclear dimers of 3d 
and 4d transition metal elements. 

While spin magnetism is routinely considered in spin 
density functional calculations, the orbital contribution 
to the magnetization is frequently neglected. This has 
two reasons: (i) in systems with high coordination like 
large clusters or bulk, the orbital moment is almost 



quenched due to chemical bonding; thus, the magne- 
tization almost exclusively originates from the spins; 
(ii) technically, orbital magnetic properties require much 
more (at least about one order of magnitude) computa- 
tional effort, as spin-orbit effects have to be considered. 
On the other hand, it is spin-orbit interaction that brings 
about one of the most important properties of magnetic 
systems: the magnetic anisotropy (MA). At low temper- 
ature, a given system resides in its magnetic ground state 
if the magnetization directs along the so-called easy axis 
of the system. The magnetic anisotropy energy (MAE) is 
the energy amount which would be needed to rearrange 
the magnetic moments of the system in the direction of 
the so-called hard axis. (Depending on the symmetry, 
several hard axes or hard planes may exist.) Evidently, 
the larger the MAE the more stable the system behaves 
against magnetic moment flips and a possible informa- 
tion stored in this way is more safe. Thus, systems with 
large MAE are interesting for technological application, 
e.g., in magnetic storage devices. 

Due to enlarged computational resources and more ad- 
vanced numerical routines, density functional methods 
are nowadays capable of reliably determine total energy 
differences well below mHartree accuracy in many cases. 
This accuracy is at least required when calculating MAE, 
which are small total energy differences. A scan through 
the existing literature reveals that the magnitude of MAE 
depends on the dimensionality of the system under con- 
sideration. As an illustrative example we refer to the 
estimates for Co systems of different dimensionalities as 
given by Gambardella et al} Therein, the MAE is esti- 
mated to increase by roughly two orders of magnitude 
when decreasing the dimensionality of the Co system. 
Amounting to 0.04 meV (per atom) for bulk hep Co, the 
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MAE increases to 0.14 meV (per atom) for a monolayer 
of Co atoms, and increases further to 2.0 meV (per atom) 
for a chain of Co atoms with averaged length of 80 atoms. 
Now, the most interesting question is, how large the MAE 
of a transition metal system of even smaller size might 
be. This leads us directly to the 3d and Ad transition 
metal dimers, being among the smallest systems showing 
anisotropy effects. We will demonstrate that potentially 
the MAE can be two further orders of magnitude higher 
than discussed by Gambardella et alr^ 

Substantial experimental work has been carried out 
on transition metal dimers. It has been summarized in 
a recent review by Lombardi and Davisi^ There is also 
an older review of Morse, ^ summarizing the experimen- 
tal data up to 1986. It has been mentioned by several 
authors that the computational investigation of dimers, 
though seemingly the most simple chemical system, is far 
from being trivial. This holds in particular for transition 
metal dimers, where frequently a number of low-lying 
states of different spin multiplicity and symmetry com- 
pete. There are several works published on the theoret- 
ical investigation of 3d and Ad transition metal dimers, 
mainly focusing on the structural properties. The 3d 
transition metal dimers have recently been investigated 
by Yanagisawa et alA using a hybrid DFT method. In 
the same year, Barden et al^ performed calculations us- 
ing six density functional or hybrid Hartree-Fock/density 
functional theory methods and an unrestricted Hartree- 
Fock method. Later, Gutsev et al^ reported on calcu- 
lations using several approximations to the density func- 
tional as implemented in the Gaussian 98 program pack- 
age, whereas Valiev et alX applied the projected aug- 
mented plane wave method to calculate structural prop- 
erties of the 3d transition metal dimers. 

Continuing with the Ad transition metal dimers, we are 
aware of only two works on their structural properties. In 
an early work, Wu^ investigated the Ad transition metal 
dimers from Y2 to Cd2 except Tc2 utilizing several ap- 
proximations to the density functional as implemented 
in the Gaussian 98 suite of programs. More recently, 
Yanagisawa et al^ reported on their investigation on Ad 
transition metal dimers from Y2 to Ag2 . They again ap- 
plied several approximations to the density functional. 
Additionally, they carried out calculations using M0ller- 
Plesset second-order perturbation method for compari- 
son. 

All the abovementioned work is restricted to the inves- 
tigation of structural properties, such as bond lengths, 
harmonic vibrational frequencies, and dissociation ener- 
gies. We are aware of only two, very recent publications 
devoted to the calculation of orbital moments and MAE 
for transition metal dimers. Strandberg et al.^^ published 
data on the MAE of C02, Ni2, and Rh2, and Fcrnandez- 
Seivane and Ferrer reported related data for Pd2, Ir2, 
Pt2, and Au2.JA Both pubHcations used local spin den- 
sity approximation (LSDA) and/or generalized gradient 
approximation (GGA) to the density functional. We are 
not aware of any MAE calculation for a complete 3d or 



Ad series of homonuclear transition metal dimers. Fur- 
ther, we are not aware of any calculation for the MAE of 
dimers beyond local approaches to the density functional. 
Both topics will be addressed in the present study. The 
presented results might become an issue of technological 
relevance once one would be able to fix the position and 
geometrical orientation of dimers without much disturb- 
ing their magnetic state. 

In the following Section [TTl we will briefly review the 
necessary theoretical background for the use of orbital 
dependent potentials within density functional theory. 
The treatment of the magnetic anisotropy in general and 
the need to introduce so-called orbital polarization cor- 
rections (OPC) is discussed in Section [1111 whereas Sec- 
tion IIVI summarizes the computational details. In Sec- 
tion IV Al the results on ground state spin multiplicity, 
bond lengths and harmonic vibrational frequencies are 
presented and discussed. Section |VB] compiles results on 
the magnetic properties (spin and orbital magnetic mo- 
ments and MAE), calculated without and with orbital 
polarization corrections. Finally, Section IVTl summarizes 
the findings of the present work and gives an outlook. 



II. ORBITAL DEPENDENT POTENTIALS IN 
DFT 

Density functional theory (DFT) has proved an enor- 
mous predictive power for real materials, although the 
model density functionals in use (LDA, GGA, B3LYP, 
. . . , just to mention a few acronyms) are partially based 
on uncontrolled model assumptions. Nevertheless, fur- 
ther progress of DFT depends on designing better model 
density functionals for chosen fields of application. (Al- 
though the existence of a universal density functional 
has been proven, there is little hope to model this prob- 
ably extremely complex functional in general.) In re- 
cent time, more and more orbital dependent function- 
als like self-interaction correction (SIC), exact exchange 
(EXX) , local spin density approximation plus local corre- 
lation (LSDA-I-C/), OPC and others came into the game. 
(Local density approximation plus dynamical mean field 
theory (LDA-I-DMFT) is a combined density-functional 
plus other-many-body-techniques approach to the self- 
energy of the single-particle Green's function and is an- 
other important development aiming mainly at electronic 
excitations.) Orbital dependence means non-locality of 
the potential operator (integral operator), and this costs 
appreciable additional computing effort. Therefore, lo- 
cal approximants of these potentials like optimized ef- 
fective potentials (OEP) are popular. Sometimes it is 
also incorrectly claimed that density functional methods 
would demand local effective potentials. True, if the uni- 
versal density functional H[n] has a unique functional 
derivative, this corresponds to a local external potential 
(w(r) — fi) = —5H/Sn{r) where fi is the chemical poten- 
tial. This results from the Hohenberg-Kohn variational 
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principle 

E[v,N\ = min{i7[7i] + (i; | n) | (1 | n) = A^}, (1) 



where 



{v\n) = 



E 



(2) 



and the spin density components vary in L^iV) of cu- 
bic summable functions on the position space V of finite 
volume.'-^ Cubic summabihty comes from the demand of 



finiteness of kinetic energy^iS The existence and unique- 
ness of the functional derivative of H[n] on hypersurfaces 
of constant particle number N or particle numbers of 
the spin channels of non-relativistic theory is the state- 
ment of the first theorem or lemma by Hohenbcrg and 
Kohn (adapted to a whole functional space for the varia- 
tion, comprising also cases of degenerate ground states). 
It is equally well known that this derivative jumps across 
those hypersurfaces. 

To perform the variation, one may introduce the gen- 
eralized Kohn-Sham ansatz 



H[72] = K[n] + L[n\, K[n\ = min {/C[(/)fe, rife] I {(t>n<j)) = n, {cf^k \ <Pi) = 4i, < < l} 



(3) 



r 



with the abbreviation orbitals 4>k{rs) and their occupation numbers nu under 

the indicated constraints, while L = H — K further de- 
n = Uss'ir) = ^cl)k{rs)nk(j)l{rs') = {<pn(f>) (4) pends on (or determines) the level of knowledge of the 
k total density functional 



and where )C[(f>kj nk] is an explicitly given functional hav- 
ing a unique minimum with respect to variation of the 



Then, the Kohn-Sham variational principle 



E[v,N] = min < K.[(j)k,nk] + L[((/>n0)] -I- y^icfiil^/rriivSik - eik)^/ri^\(j)k) < < 1, 'S^ nk = n\ 



(5) 



r 



is obtained. The eik are Lagrange parameters for the 
orthonormality constraints. If as usually K, includes an 
kinetic energy term of orbitals, then the minimal density 
n from ([5]) is automatically in L^{V). With the writing 



Uk Scj^l 



fk[<Pi,ni] 
rik 



Pk), 



SL 
Sn 



(6) 



the variation leads to the generalized Kohn-Sham equa- 
tions 



((/>;|fcfc VL + v\(j)k) = y/ni/ukeik 



(7) 



and to the ordinary aufbau principle for the occupation 
numbers Uk ■ In general they may be solved together with 
the orthonormality constraints by a unitary triangulation 
of the matrix of t he left hand side expressed in a suitable 
order of an orthonormal orbital basis. If as in many cases 
kk is a Hermitian operator, a diagonalization is possible. 

In the ordinary Kohn-Sham approach kk consists of 
the orbital kinetic energy operator which is the same dif- 
ferential operator for all orbitals. (The classical Coulomb 
potential of the density n can likewise be considered as 



an orbital expression or as part of wl, here we adopt 
the latter alternative.) Nevertheless, discontinuities of 
the functional derivative of K[n] when varying the par- 
ticle number are incorporated in this approach. They 
yield a piecewise linear behavior of this contribution to 
the total energy as a function of particle number and 
an appreciable part of band gap energies. Improvement 
of ionization potentials, electron affinities and band gaps 
can be expected from further orbital related terms intro- 
duced in kk as is the case in the SIC, LSDA+U, EXX, 
OPC and other approaches. In particular in magnetism, 
the non-relativistic total energy may be piecewise linear 
even with the total particle number fixed, if it changes 
between the two spin channels as analyzed in Ref. Il4l . 
Note that every piecewise linear behavior of the total en- 
ergy functional causes a discontinuity in the derivative of 
its Legendre transform which is the density functional. 
Therefore, there may be good reasons not to remove the 
orbital dependence of potentials from model functionals. 

The numerical results in the present paper are obtained 
with the LSDA and the LSDA+OPC functionals. Since 
we are aiming at magnetic anisotropy energy which is due 
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to spin-orbit coupling, the Dirac-Kohn-Sham equations 
are solved. Then, in the LSDA case, kk reduces to the 
Hermitian relativistic kinetic energy operator t and the 
orbitals arc four-component bispinor orbitals. In the case 
of LSDA+OPC one has 

kk=i+voPC, (8) 
with the Hermitian non-local potential 

^opc= I'^iC-'i^'l- (9) 

The multi- index v stands for v = (i?, p, k, ii), where R is 
the atomic position, p is the principal atomic site orbital 
quantum number, k is the spin-orbit quantum number 
linked to the orbital angular momentum quantum num- 
ber by 1{k) = K for «; > and 1{k) — — k — 1 for k < 
in the absence of external magnetic field (in the presence 
of a magnetic field this connection still holds to a very 
good approximation for the large bispinor component), 
and /i is the z-component of the total angular momentum 
with the z-axis chosen in direction of the site magnetic 
polarizationii^ The non-orthogonal four-component basis 
orbitals for the algcbraization of the Dirac-Kohn-Sham 
equation used in the full-potential local-orbital (FPLO) 
packagei^iii are denoted by Iv), and \i^} are the states 
of the contragradiant basis, {i'Ij^') = 6^,^'. The matrix 
elements of the used OPC versioni^ are 

a.' =S,(«)M,(i.|L,|i/') , (10) 

where Lz is the operator of the z-component of the orbital 
angular momentum, 

= J2MMp>}{'^\LzW')W\pM (11) 

k 

is the z-component of the obtained orbital momentum in 
the considered atomic spin-subshell (P^- being projection 
operators to the spin channel) , and Bi is a Racah coeffi- 
cient {B2 = (9F2 - 5F'*)/441 in terms of Slater integrals 
with the related 3d or Ad radial functions). The v and i'' 
sums in ^ run over the 'magnetic' basis orbitals which 
are the 3d or 4d basis orbitals in our case being quite lo- 
cal and well representing the atomic d shells of the atoms 
in the dimer (or in other cases in a solid). 



III. MAGNETIC ANISOTROPY 

Magnetic anisotropy is the dependence of measurable 
quantities on the direction of the applied magnetic field, 
B. In a more narrow sense, but much more convenient 
for theoretical considerations, one can consider MA as 
the dependence of free or inner energy on the direction 
of magnetization, M. 



In extended systems, both viewpoints are commonly 
related by splitting^- the generating function of free en- 
ergy, $(r,B,M): 

$ = $o + ^MA(e)-|M|eB. (12) 

Here, the dependence on the direction vector of magne- 
tization, e — M/\M\, and on applied magnetic field has 
been made explicit, but the first two terms also depend 
on \M\. They further depend on temperature T, pres- 
sure, etc. Note, that the subdivision into the isotropic 
part $0 and the anisotropic part Ema is not unique and 
that the free energy is the minimum of $ with respect to 
M. 

The suggested notation makes sense if the field-induced 
change of ($0 + £'MA(e)) via variation of \M\ is much 
weaker than the change of i?MA(e) via variation of e. 
This is usually the case, e.g., for ferromagnets well below 
the magnetic ordering temperature in moderate applied 
fields^ 

For the present purpose, we will adopt a somewhat 
different notation. Considering finite systems, thermo- 
dynamic features can be dropped. The focus is now on 
the dependence of ground state energy, Eq, on the di- 
rection of magnetization, to be evaluated by constrained 
DPT. Define 

MAE(e) = So(e)-So(0,0,l), (13) 

where (0,0,1) denotes the direction parallel to some 
appropriately chosen axis, e.g., the dimer axis. This 
definition does not rely on a weak dependence of the 
anisotropic part of the energy on the magnitude of the 
total moment, since the relation with previous notation is 
Eo{e) = i?MA(e, |iW|(e),r 0) for appropriate choice 
of $0. As we will see in the results, the ground state 
value of the total magnetic moment, 

j d\\M\ = \{2S + L)\p^= Ps + m (14) 

may indeed substantially depend on e. This direction is 
imposed on the system by constraining the xc field paral- 
lel to chosen e. As only coUinear spin- and magnetization 
densities are considered, 

lJis + ^J■l = e{2Sz+ Lz)p^ ■ (15) 

Note, that pi — p-Q M^- relates the notation in terms 
of magnetic quantum numbers, Eq. (jlip . with that in 
terms of orbital moments. 

Dimers are the smallest chemical objects that show 
MAfi^ While Eo{e) of atoms without external magnetic 
field is isotropic and obtains axial symmetry only by ap- 
plication of a field, dimers possess axial symmetry at 
B ~ 0. One could argue that free dimers would always 
orient themselves in the external magnetic field such that 
the lowest possible energy is achieved. This statement is, 
however, valid for any finite system. For example, macro- 
scopic samples are frequently positioned free to rotate in 
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magnetic measurements in order to ensure orientation of 
M along B.'^'^ What we aim at, however, is to investigate 
dimers with a rotation axis fixed in space by an unspeci- 
fied interaction, e.g., as a model for dimers on a substrate. 
This model interaction should be strong enough to with- 
stand the torque due to the magnetic field. On the other 
hand, it should not essentially influence the electronic 
structure of the dimer and, thus, the MAE defined in 
Eq. (O. 

One challenge to both computational and applied 
chemistry would be to find an embedding medium fuUfill- 
ing the two opposing requirements in order to verify the 
enormous MAE values predicted by our calculations and 
two preceding papers J^i^i If such a medium was found, 
the MAE values would represent barrier heights against 
thermal and quantum mechanical fluctuations of the 
super-paramagnetic dimer moments. A 10-year stability 
of super-paramagnetic bits would require MAE/fcT > 40, 
k being the Boltzmann constant.^^ Our results will show 
that such a stability could be achieved for T w 60 K. 
Thus, the search for dimers (or other very small clusters) 
with room-temperature stability should not be consid- 
ered an unrealistic goal. 

In a non-relativistic theory, spin and real space are in- 
dependent of each other and spin magnetization can be 
rotated by the external field without change of energy. 
Thus, the consideration of spin-orbit coupling (and other 
relativistic effects on the level of the Dirac-Kohn-Sham 
formalism) is needed for the description of MAE. Here, 
a four-component relativistic approach is usedfi^ Be- 
yond these kinematic relativistic effects, magnetic dipole- 
dipole interaction contributes to MAE. This interaction 
is the source of so-called shape anisotropy in macroscopic 
theory. It is responsible for the fact that iron magnets 
should have special shape, e.g., that of horse shoes or long 
bars, in order to perform their operation. In formal the- 
ory, this interaction is described as Hartree-like approx- 
imation to the Breit interaction of relativistic quantum 
electrodynamicsj22, Here, we neglect the shape anisotropy 
as it is small (about 1 meV) compared with the spin-orbit 
generated anisotropy in most of the considered dimers. 

Orbital dependent exchange effects to the MAE of 
dimers are more important in some cases and have been 
neglected in the previous approaches f^Siii Such effects are 
not included in the standard DFT approximations like 
LSDA or GGA, which are based on homogeneous elec- 
tron gas theory. For example, LSDA does not account for 
the second Hund rule (maximum orbital quantum num- 
ber in the ground state) in atoms. ^'^ Thus, it is likely 
that also the ground state of dimers is not accurately 
enough described for the present purpose in any known 
local approximation. For this reason, we have decided to 
add orbital polarization corrections (OPC) to LSDA, as 
described in the previous section, in the present calcula- 
tions. Such corrections were first suggested by Brooks, 
Eriksson, and Johansson] ^^'^^ The specific implementa- 
tion used is taken from Ref. [l^. Though this is just 
one out of several similar approaches it will provide an 



impression how strong the influence of orbital exchange 
terms on such a sensitive quantity like MAE can be. As a 
matter of (scarce) experience, LSDA gives a lower bound 
on the MAE values while LSDA-f OPC gives an upper 
bound, at least for systems with low coordination. An 
example where both experimental and theoretical data 
are available is Co atoms and small Co clusters on the 
surface of Pt. Here, the strength of OPC had to be down- 
scaled by 50% in order to meet the experimental values of 
MAE.^^ Below, we will present and discuss both LSDA 
and LSDA-I-OPC data. This will allow the reader to get 
an idea of the order of magnitude that can be expected 
in reality. 

While DFT approaches to (orbital) magnetic proper- 
ties of extended systems are state-of-the-art, considerable 
less experience is available for small clusters. We thus feel 
inclined to provide the following arguments to justify our 
approach: 

• Spin magnetism and spin-orbit coupling in the 
ground state of single atoms and of weakly corre- 
lated solids are very well described in LSDA.— 

• Orbital magnetism in the ground state of sin- 
gle atoms is by construction well described in 
LSDA-I-OPC; it works reasonably well in metallic 
magnets 

• There are cases (e.g., lanthanides) where the 
ground state configuration and ionization poten- 
tials of atoms are not well described in LSDAi^ 
This problem is not present in bulk metals. 

• As a summary of the three points above, one can 
expect relativistic LSDA(-I-OPC) to work reason- 
ably well for dimers, if one can assure the correct 
ground state configuration and ground state spin 
multiplicity. Therefore, we compare the results of 
scalar-relativistic calculations with available results 
of other calculations and with experimental data. 

• Last but not least, we are aware of only two LSDA 
and/or GGA calculationsiSiii to study orbital mag- 
netic properties, in particular MAE, of dimers. No 
LSDA-I-OPC or more involved quantum chemical 
calculations on this topic have been reported so far. 
Thus, the presented results are the first that will 
not only allow to estimate a lower bound (LSDA 
and GGA) but also an upper bound to the mag- 
nitude of MAE and trends among the considered 
series of transition metal dimers. 

Before we proceed to more computational details, a fi- 
nal remark on the angle dependence of MAE is in place. 
In dimers, MAE, Eq. due to axial symmetry only 

depends on the angle O between the dimer axis and the 
magnetization. Further, extrema of MAE(O) are ex- 
pected at O = and 9 = 7r/2, for symmetry reasons. 
Other extrema cannot be excluded but are unlikely. They 
could occur, e.g., due to multiple level crossing upon vari- 
ation of Q. The present investigation is constraint to the 
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evaluation of MAE(7r/2), MAE(O) being zero by defini- 
tion. A more detailed investigation should scan MAE(8) 
at some finer mesh. 



IV. COMPUTATIONAL DETAILS 

We are aiming at the MAE which is due to spin-orbit 
coupling and hence of the order of the fine structure split- 
ting of the molecular terms. Apart from exceptionally 
possible cases, the distance in energy between molecular 
terms at equilibrium bond lengths (typically several 100 
meV) is larger than the fine structure splitting, so one 
may concentrate on the MAE related to the molecular 
ground state term, where for the latter spin-orbit split- 
ting has been neglected. After having identified this term 
and its structural properties, spin-orbit coupling will be 
switched on to study MAE. 

Without spin-orbit coupling, a molecular term ^'^+^Ap 
is characterized by the spin multiplicity 25* -1-1 of the total 
molecular spin S, by the absolute value of the eigenvalue 
M of the projection of the total orbital momentum on the 
molecular axis (A H, A, . . . for \M\ = 0, 1, 2, . . .), 

and by the parity p = g^u with respect to inversion 
in electron configuration space. (E^ are M = states 
which are even/odd with respect to reflection on a plane 
containing the molecular axis; an M 7^ state is trans- 
formed into a —M state of the same term by such a 
reflection.) Besides the spin multiplicity, there is an or- 
bital multiplicity. Within the scalar-relativistic approxi- 
mation (spin-orbit coupling neglected), spin and orbital 
parts of the state are only related via fermionic antisym- 
metry with respect to particle exchange, and the spin 
density matrix with respect to a chosen z-axis of spin may 
be well defined in the states of each molecular term since 
the Hamiltonian is spin independent. It may be chosen 
either even or odd with respect to an inversion through 
the center of the molecular geometry (sometimes also 
called symmetric and symmetry-broken ground state). 
In order to obtain single-determinant Kohn-Sham states 
representing a given term, coUinear fixed spin moment 
(FSM) calculationsSS. are performed with initially even 
spin density and predefined Sz — S. For S = 0, alterna- 
tively odd spin density is enforced by initial conditions 
for the self-consistency iterations. The latter case corre- 
sponds in particular to the ground state for Cr2, Mn2, 
and M02. Calculations for all S*- values are performed, 
and the ground state term is determined as that of low- 
est total energy. 

These calculations are performed with the scalar- 
relativistic version of the full-potential local-orbital code 
FPLQj ^^i^^ Its newly implemented feature for calculat- 
ing the electronic structure of clusters or molecules (up 
to several 100 atoms) is used. The valence basis com- 
prised 3s3p3d 4s4p4(i 5s states in the case of 3d transi- 
tion metal dimers, while a basis with the respective next 
higher principle quantum number was used for the 4d 
case. Within the LSDA, the exchange-correlation func- 



tional by Perdew and Wang (PW92)^° is applied. A fur- 
ther trick within the FSM to improve the convergence be- 
havior has been the variation of the temperature broad- 
ening parameter for the Brillouin zone integration (from 
500 K stepwise down to 100 K), and varying the nu- 
merical procedures for the self-consistency loop within 
the calculations. The obtained ground state multiplici- 
ties 25* -|- 1 for the 3d and Ad transition metal dimers are 
given in the second column of Tabs. HI and [Til respectively. 

After having identified the ground state spin multi- 
plicity and symmetry for each of the 3d and Ad transi- 
tion metal dimers within our LSDA approach in scalar- 
relativistic approximation, we have obtained the bond 
lengths and harmonic vibrational frequencies by analyz- 
ing the total energy curves belonging to the dimer ground 
states. The necessary curve fitting has been done utiliz- 
ing cubic splines. Exemplarily, the total energy curves 
for various spin multiplicities 25* -I- 1 are given for the 
C02 dimer in Fig. [TJ 

Being mainly interested in the orbital magnetic prop- 
erties of the 3d and 4d transition metal dimers, we have 
then performed fuU-relativistic FPLO calculation a^^'^^ at 
the dimer ground state bond lengths and spin multiplic- 
ities obtained by previous calculations. It is known from 
earlier investigations on transition metal bulk systems, 
see discussion in Sec. IIIIl that treating the orbital mag- 
netism in LSDA or GGA leads to severe disagreement 
with experimental observations. Therefore, additionally 
orbital polarization corrections have been introduced to 
calculate orbital moment related properties. Within this 
work all the 3d and id transition metal dimers have been 
treated within the standard fuU-relativistic LSDA ap- 
proach. Additionally, we performed calculations taking 
into account orbital polarization corrections, as described 
m Section [TTl 



V. RESULTS 

In this section, collecting the results of the calculations, 
we first give the structural properties of the dimers for the 
ground state spin multiplicities 25+1 obtained within 
the scalar-relativistic LSDA approach. In a second part, 
all results on the magnetic properties obtained by fuU- 
relativistic calculations are summarized. 



A. Spin multiplicities and structural properties 

The bond lengths and harmonic vibrational frequen- 
cies for the 3d and Ad transition metal dimers obtained 
within the scalar-relativistic LSDA approach are given 
for the ground state spin multiplicities in Tables HI and HIl 
respectively, and shall now be discussed in detail. The 
tables are arranged such that the first line for each dimer 
contains our ground state data and other data which re- 
fer to the same ground state. If other calculations or 
experiment yield another ground state, additional lines 
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are added and our results for this state (being an ex- 
cited state according to our calculation) are included for 
comparison. 

Starting with the 3d transition metal dimers, our ap- 
proach predicts the ground state of the early transition 
metal dimers Sc2, Ti2, and V2 to have spin multiplicity 
5,3, and 3, respectively, in agreement with other calcula- 
tions and experimental observations. The bond lengths, 
given in Fig. [2] (a), agree reasonably well (within 2%) 
with other LSDA calculations performed by Barden et 
al.^ An even better agreement (within 1%) is obtained, 
probably by coincidence, with the B3LYP calculations 
by the same authorsr A similar good agreement is found 
between our calculations and experiment for Ti2 and V2 
(note, that the value for Sc2 given in column 'Exp' is 
an estimation). The calculated vibrational frequencies 
agree well with each other and deviate in a systematic 
way from the experimental data, see Fig. [2] (b). 

Continuing further, the Cr2 dimer shows competition 
between low-lying even and odd zero-spin states. Our cal- 
culation yields an odd ground state, in agreement with 
the LSDA calculation by Valiev et alJ This is at vari- 
ance with the results by Barden et aZ.— who report an 
even ground state obtained both in LSDA and in B3LYP. 
We show structural data for both states in Table [D In 
our calculation, the even state (with zero spin density) 
lies 127 meV above the symmetry-broken ground state. 
The respective bond lengths agree within about 2% as for 
the previous examples. Remarkably, the vibrational fre- 
quencies are considerably different for both states, the 
state with a non-zero but odd spin density providing 
a much softer harmonic potential than the zero spin- 
density state. As the experimental frequency agrees with 
the calculated frequencies belonging to the odd state, we 
assume that this is the true ground state and present the 
experimental data in this line. The Cr2 dimer is excluded 
from the further investigation, since its zero-spin ground 
state is in any case hardly susceptible to moderate ap- 
plied external magnetic fields. 

The Mn2 dimer is difficult to treat within LSDA. This 
is due to the high promotion energy from the atomic 
(3(i)^(4s)^ configuration to the (3c?)^(4s)^ configuration 
preventing strong hybridization and making van der 
Waals interactions important Further, several ground 
states with different spin multiplicities between 1 and 11 
compete. We find a ground state with the latter value in 
agreement with DFT calculations performed by Gutsev 
et ali^ However, the experimental ground state is a sin- 
glet state, probably with odd non-zero spin dcnsityi^i^ 
The Mn2 dimer is therefore excluded from further con- 
sideration. 

For the remaining 3d transition metal dimers we ob- 
tain again good agreement with other calculations and 
fair agreement with experiment. No experimental data 
are known for the ground state spin multiplicities of Fe2 
and C02. Here, our results agree very well with other 
calculations. For Fe2, also the experimental bond length 
compares well with the calculations, but for both Fe2 and 



C02 the vibrational frequencies are lower in the experi- 
ment than in the DFT calculations. The same holds for 
Ni2, where the experimental ground state is reported to 
be a mixture of a triplet and a singlet state due to spin- 
orbit coupling.^ Our calculations yield a triplet ground 
state, in agreement with the LSDA calculation by Valiev 
et al.^ but in disagreement with the calculation by Bar- 
den et al^ who find a singlet state with even spin density. 
We also provide data for this state which is found at 414 
meV in our calculation. The full-relativistic calculations 
will be carried out for the triplet state of Ni2. They 
should be taken with caution in view of the mentioned 
spin-orbit mixing between singlet and triplet states. 

The copper dimer has a singlet ground state with zero 
spin density and is not considered further. Its structural 
properties obtained in the calculations agree reasonably 
well with the experimental data. 

Continuing with the Ad transition metal dimers, their 
bond lengths and harmonic vibrational frequencies are 
given for the respective ground state multiplicities in 
Tab. [ni Quite similar to the id transition metal dimers 
discussed above, the results obtained for the early Ad 
transition metal dimers Y2, Zr2, and Nb2 are in agree- 
ment with other theoretical calculations and experiment. 
The correct ground state spin multiplicities of 5, 3, and 
3, respectively, are obtained for these systems in our cal- 
culations. The bond lengths and harmonic vibrational 
frequencies agree with other reported values to within a 
few percent, as can be seen in Fig. [3] (a) and (b). The 
FPLO results for the Ad transition metal dimers are in 
a better agreement with the calculations performed by 
Yanagisawa et alA than with the calculations performed 
by Wu.? 

Proceeding further, the Cr column transition metal 
dimer M02 has been reported to show a broken-symmetry 
(odd zero-spin) ground state in LSDA.'^^ Our calculation 
yields a zero spin density ground state, but the odd zero- 
spin ground state lies only a few meV above. We exclude 
M02 from further consideration. Experimental and theo- 
retical structural data show reasonable agreement, with 
a somewhat smaller vibrational frequency in experiment 
compared with the calculations. One is tempted to take 
this difference as hint to an odd zero-spin ground state 
in experiment. On the other hand, the vibrational fre- 
quency calculated for such a state, 520 ± 50 cm"-'^,^^ is 
not clearly different from the related frequencies calcu- 
lated for even zero-spin ground state, at variance with 
the situation in Cr2. 

For the Ru2 dimer, we obtain a ground state spin mul- 
tiplicity of 5, in disagreement with other calculations 
which yield a ground state multiplicity of Our cal- 
culated state with 25 -I- 1 = 7 lies approximately 400 
meV above the ground state. Its structural data agree 
well with related literature values and also with the ex- 
perimental vibrational frequency. The latter supports a 
25 + 1 = 7 ground state. While we could stabilize this 
state in the scalar-relativistic calculations, we failed to 
achieve this in the full-relativistic approach. This is the 
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reason why no data for Ru2 are provided in the next 
section. 

For the remaining Ad transition metal dimers Rh2, Pd2, 
and Ag2 we can adopt the discussion of the late 3d tran- 
sition metal dimers. The agreement with other calcu- 
lations and experimental results is reasonable, and the 
agreement with calculations performed by Yanagisawa et 
ai^ is better than with results by Wu.^ The non- magnetic 
Ag2 dimer is excluded from further consideration. 

In general, we find quite good agreement for the 3d 
and 4c? transition metal dimers concerning the ground 
state spin multiplicity, their bond lengths and harmonic 
vibrational frequencies. Ground states with wrong spin 
multiplicity are found for Mn2 and Ru2 . Comparing our 
results with the most recent work by Strandberg et al.^ 
who considered C02, Fe2, Ni2, Rh2, Ru2, and Pd2, we 
find the same ground state multiplicities for all systems 
including Ru2. 



B. Orbital magnetic properties 

After having identified the spin multiplicities of the 
dimer ground states within a scalar-relativistic treat- 
ment, we now discuss the results on the orbital magnetic 
properties obtained by the fuU-relativistic approach to 
the dimer ground states. The spin (/i^) and orbital (/x;) 
moments of the 3d and 4d transition metal dimers are 
given in Tables IIIII and IIVI respectively, for magnetiza- 
tions directing perpendicular to and along the dimer axis. 
For each dimer the results of the fuU-relativistic LSDA 
calculations are given in the first row, whereas in the sec- 
ond row results are compiled obtained by fuU-relativistic 
LSDA-I-OPC calculations as discussed in Sec. IIIII Gen- 
erally, in small systems like dimers orbital moments are 
less quenched compared to systems with larger coordina- 
tion numbers, and second Hund's rule couplings are more 
important. 

The orbital moments are relatively small for the first 
and third transition metal dimers of the 3d and 4d rows. 
This is in agreement with the experimentally observed^ 
zero orbital angular momentum ground states. (The val- 
ues of the ground state orbital angular momenta com- 
piled by Lombardi et al.^ are given in parentheses behind 
the easy-axis calculated orbital moments. Tables lllll and 
IIV[) . Orbital polarization corrections enhance the orbital 
moments up to a factor of two, but the MAE is rela- 
tively small for both LSDA and LSDA+OPC approaches, 
as spin-orbit coupling has little influence on the L = 
states. 

The situation is different for the Ti2 and Zr2 dimers 
with experimentally observed L = 2 ground states. Here, 
we obtain significant orbital moments and anisotropy en- 
ergies. As a grain of salt, the evaluated orbital moment 
of the Ti2 dimer is smaller than expected. This might 
be a problem of LSDA. It is only partly cured by the 
application of OPC, which changes the sign of the or- 
bital moment along the dimer axis of Ti2. This is related 



with a reduction of the MAE. The opposite holds for Zr2, 
where OPC enhances the MAE to a considerable value 
of 0.08 eV. In this case, the calculated orbital moment is 
close to the experimental one. It should be noted that 
these experimental values exist only for the easy axis of 
magnetization, i.e., for the ground state orientation. The 
same holds for Fe2 and C02, where the orbital moments 
along the easy axis are very close to 2 reported in 
literature (calculations), but are much smaller along the 
hard axis. For Ni2 and, partly, for Pd2 the situation is 
less clear. This might be due to a stronger spin-orbit 
mixing in these systems and possible level crossings that 
deserve further investigation. 

Very interesting results are obtained for the magnetic 
anisotropy energy of the late transition metal dimers be- 
ginning with Fe2 and Tc2 for the 3d and 4d row, respec- 
tively. For the C02 and Rh2 dimers we obtain an MAE of 
-1-50 meV (LSDA) and +104 meV (LSDA), respectively, 
in agreement with -f 60 meV and -t-90 meV reported by 
Strandberg et alt^ Due to OP corrections these values 
are enhanced to -1-188 meV and -1-196 meV, respectively. 
Another large MAE value is found for Fe2, almost com- 
peting in size with the two previously mentioned dimers, 
C02 and Rh2. The record value of 223 meV obtained for 
Ni2 by LSDA-f OPC should be taken with caution, since 
Ni2 has a multiconfiguration ground state. 

Intriguingly, all large values of MAE in the order of 0.1 
eV are related with uniaxial anisotropy and anisotropy of 
the orbital moment magnitude. This can be understood 
from the interplay of bonding and spin-orbit splitting. If 
the magnetization is parallel to the dimer axis, orbitals 
with zero magnetic quantum number form a bonds. The 
states with highest magnetic quantum numbers, \mi \ = 2, 
form 5 bonds instead. Spin-orbit interaction can gain 
considerable energy by splitting these states. On the 
other hand, if the magnetization and, thus, the quan- 
tization axis, is perpendicular to the dimer axis, the a 
bonds are formed by combined states with to/ — ±2. In 
this case, there is no effect of spin-orbit interaction on 
these states and the energy is higher. Similar arguments 
have been provided by Strandberg et al<^ 



VI. SUMMARY AND OUTLOOK 

We have reported a systematic density functional 
study of ground state properties, in particular of the 
magnetic anisotropy energy, for the homonuclear 3d and 
4d transition metal dimers. The notation of magnetic 
anisotropy, central to this paper but less common in 
chemistry, has been introduced in some detail. Before, 
orbital dependent potentials in density functional theory 
have been discussed, as a reasonable description of orbital 
magnetism including magnetic anisotropy requires cor- 
rections to the most common LSDA or GGA approaches 
within DFT. The presented results have been obtained by 
scalar- and fuU-relativistic LSDA calculations, the latter 
without and with orbital polarization corrections applied. 
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To ensure consistency of the calculations with previous 
work, ground state spin multiplicities, bond lengths, and 
vibrational frequencies have been compared with avail- 
able theoretical and experimental data from literature. 
Good agreement between the present data and results 
from literature has been found for most cases. Spin and 
orbital moments as well as MAE have been evaluated for 
all dimers with 25* + 1 > 1 in the ground state, except 
Ru2, where our approach did not find the correct ground 
state spin multiplicity. 

The main finding of the present work is an excep- 
tionally large magnetic anisotropy energy of about 0.2 
eV/f.u. for Fe2, C02, Ni2, and Rh2, obtained by means 
of LSDA-hOPC. Values of 0.06 eV/f.u. and 0.09 eV/f.u. 
have recently been reported in literature-'^'^ for C02 and 
Rh2, respectively, obtained by means of LSD A calcu- 
lations. We confirm these LSDA values but note that 
LSDA is known to provide, in many cases, only a lower 
bound on the strength of orbital magnetic properties in- 
cluding MAE, while LSDA-I-OPC frequently yields re- 
sults closer to experiment. '^^^ Another set of LSDA and 
GGA calculations on dimer MAE was reported very re- 
cently, where large numbers have been found for the 
MAE of Pt2 (0.035 eV with LSDA and 0.1 eV with 
GGA) and of Ira (0.06 eV with LSDA).ii The quoted 
paper also reports an LSDA value, -0.005 eV, for the 
MAE of Pd2, in agreement with our related value. Fi- 
nally, a small negative MAE is found for Ni2 in Ref. 
fiol . in disagreement with our positive value of 0.011 eV. 
We attribute this difference to the mentioned problem 
to describe the ground state of Nia with LSDA-like ap- 



proaches. A fuU-relativistic many-electron treatment of 
such systems would be desirable but seems hard to per- 
form at present. 

All these data raise expectations about the potential 
applicability of magnetic dimers for information storage. 
In particular, a barrier height of 0.2 eV as obtained in 
the present study would translate into 10 years stability 
of superparamagnetic bits at a temperature of 60 K. It 
is therefore tempting to search for other homo- or het- 
eronuclear dimers with yet higher MAE. A further nec- 
essary and perhaps not unrealistic step would be to find 
a matrix that can fix the dimers in space without deteri- 
orating their magnetic properties. Finally and probably 
most difficult, a completely new read/write technology 
would have to be invented for any practical application 
of such small magnetic entities. 

Besides these dreams, exceptional magnetic properties 
have always been a playground for fundamental research. 
Magnetic dimers are potential model systems in this field 
due to their unique combination of structural simplicity 
and large magnetic anisotropy. 
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FIG. 1: Total energy curves of the C02 dimer obtained within 
scalar-relativistic approximation for various spin muhiplici- 
ties 2S + 1. 
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Sc^ Ti^ Cr^ Mn^ Fe^ Co^ Ni^ Cu, 

FIG. 2: Structural properties of the 3d transition metal 
dimers: (a) equilibrium bond length do [A], and (b) harmonic 
vibrational frequency cje [cm~^]. Note, that some estimated 
values are included for the experimental bond lengths and 
some data refer to ground states of different spin multiplici- 
ties. More details are given in the related Table |T] 
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FIG. 3: Structural properties of the 4d transition metal 
dimers: (a) equilibrium bond length do [A], (b) harmonic 
vibrational frequency tOe [cm~^]. Note, that some estimated 
values are included for the experimental bond lengths and 
some data refer to ground states of different spin multiplici- 
ties. More details are given in the related Table HTl 
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TABLE L Spin multiplicity 25' + 1, bond length do [A], and harmonic vibrational frequency ujf, [cm~^] for the homonuclear 
3d transition metal dimers. Results from this work (abbreviated TW) are compared with experiment and with other DFT 
calculations from literature. Zero-spin states with even (zero) spin density are denoted by (e) and those with odd (non-zero) 
spin density by (o). 



Dimer 2S + 1 do [A] uOe [cm ] 





TW 


Exp" 


TW 


lsda' 


B3LYP'' 


Exp" 


TW 


lsda' 


B3LYP" 


Exp" 


SC2 


5 


5 


2.598 


2.569 


2.595 


(2.46)'' 


270 


263 


263 


240 


Tia 


3 


3 


1.904 


1.861 


1.910 


1.942 


488 


533 


498 


408 


Va 


3 


3 


1.765 


1.741 


1.747 


1.77 


679 


687 


687 


537 


Cra 


l(o) 


1 


1.700 


1.66= 




1.679 


448 


557" 




481 




1(e) 




1.626 


1.607 


1.605 




799 


821 


837 




Mna 


11 




2.524 




2.663 




224 




195 






l(o) 


1 


2.429 






3.4 


144 






76.4 


Fea 


7 




1.964 


1.964 


1.991 


2.02 


460 


442 


435 


300 


Coa 


5 




1.946 


1.937 




(2.16)'' 


424 


426 




297 


Nia 


3 


3« 


2.051 


2.054" 




2.154 


367 


328" 




259 




1(e) 




2.054 


2.058 


2.296 




360 


347 


227 




Cua 


1(e) 


1 


2.161 


2.179 


2.281 


2.219 


298 


289 


240 


266 



"Experimental results taken from Lombardi and Davis— 
'LSDA calculations by Barden et al^ 
'^B3LYP calculations by Barden ei al^ 
''Estimated in Ref . H using Pauling's relation. 
"Taken from calculations by Valiev et al^ 

■^Suggested from density functional calculations quoted in Ref. [^. 
*Due to strong spin-orbit coupling, the ground state is a mixture 
of states with spin multiplicities 3 and 1 — 



TABLE II: Spin multiplicity 25* -f- 1, bond length do [A], and harmonic vibrational frequency [cm~^] for the homonuclear 
4ci transition metal dimers. Results from this work (abbreviated TW) are compared with experiment and with other DFT 
calculations from literature. Zero-spin states with even (zero) spin density are denoted by (e) and those with odd (non-zero) 
spin density by (o). 



Dimer 25* -H do [A] [cm ] 





TW 


Exp" 


TW 


LSDA' 


B3LYP" 


Exp" 


TW 


LSDA*" 


B3LYP" 


Exp' 


Ya 


5 


5 


2.911 


2.716 


2.907 


(2.62)^* 


188 


207 


184 


184 


Zra 


3 


3 


2.273 


2.439 


2.248 


2.241 


309 


298 


335 


306 


Nba 


3 


3 


2.115 


2.245 


2.064 


2.078 


465 


363 


450 


425 


Moa 


1(e) 


1 


1.972 


1.977 


1.976" 


1.929 


550 


546 


548" 


477 


Tea 


3 


{3}" 


1.994 




1.928 


(2.10)'' 


526 




558 




Rua 


5 




2.114 








425 












{^y 


2.238 


2.236 


2.274 


(2.17)'' 


348 


327 


306 


347 


Rha 


5 


{^y 


2.216 


2.213 


2.203 


(2.29)'' 


351 


352 


311 


284 


Pda 


3 


3 


2.437 


2.934 


2.503 


(2.48)'' 


234 


140 


206 


210 


Aga 


1(e) 


1 


2.506 


2.505 


2.596 


2.530 


207 


212 


176 


192 



"Experimental results taken from Lombardi and Davis™ 

'LSDA calculations by Wu.? 

"B3LYP calculations by Yanagisawa et aim- 

''Estimated in Ref. using Pauling's relation. 

"Taken from calculations by Wuc§^ 

^Suggested from quantum chemical calculations quoted in Ref. [^. 
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TABLE III: Spin /is and orbital /i; moments of the 3d transition metal dimers obtained by fuU-relativistic calculations with 
magnetic moment orientation perpendicular and parallel to the dimer axis, respectively. The different values obtained by LSDA 
calculations and by LSDA+OPC calculations as discussed in Sec. lIIII are given in the first and second row for each of the dimers, 
respectively. MAE(7r/2) is the magnetic anisotropy energy whose sign is defined in Eq. (|13p . MAE(7r/2) = -Eo(-L) — i-odl), 
i.e., positive values correspond to an easy axis parallel to the dimer axis. All values are given per dimer. Moments deviating 
less than 0.01 fj,B from an integer value are rounded to the next integer. The values in parentheses given behind the easy-axis 
orbital moments are experimental {theoretical} values of the corresponding ground state orbital angular momentum given in 
Table 2 of Ref. 2. 



Dimer 


/is 


m 


/is 




MAE(7r/2) 




[/is] 


[/is] 


[fJ.B] 


[fJ-s] 


[meV] 


SC2 


4. 


-0.040 (0) 


4. 


0. 


- 




4. 


-0.062 (0) 


4. 


0. 


- 


Ti2 


2. 


-0.016 


2. 


-0.264 (2) 


+ 36 




2. 


-0.026 


2. 


0.434 (2) 


+ 33 


V2 


2. 


-0.016 (0) 


2. 


0. 


- 1 




2. 


-0.030 (0) 


2. 


0. 


- 1 


Fe2 


6. 


0.186 


6. 


1.890 ({2}) 


+ 32 




6. 


0.402 


6. 


1.908 ({2}) 


+ 150 


C02 


4.144 


0.334 


4.088 


2. ({2}) 


+ 50 




4.142 


1.196 


4. 


2.096 ({2}) 


+ 188 


Ni2 


1.986 


0.454 


1.988 


0.878 


+ 11 




1.990 


1.468 


2. 


3.120 


+ 223 



TABLE IV: Spin /is and orbital moments of the 4d transition metal dimers obtained by fuU-relativistic calculations with 
magnetic moment orientation perpendicular and parallel to the dimer axis, respectively. The different values obtained by LSDA 
calculations and by LSDA-(-OPC calculations as discussed in Sec. lIIII are given in the first and second row for each of the dimers, 
respectively. MAE(7r/2) is the magnetic anisotropy energy whose sign is defined in Eq. (|13|l . MAE(7r/2) = -Eo(-L) — -Eo(||), 
i.e., positive values correspond to an easy axis parallel to the dimer axis. All values are given per dimer. Moments deviating 
less than 0.01 fiB from an integer value are rounded to the next integer. The values in parentheses given behind the easy-axis 
orbital moments are experimental {theoretical} values of the corresponding ground state orbital angular momentum given in 
Table 2 of Ref. '5. 



Dimer 


/is 




/is 


l-H 


MAE(7r/2) 




[/is] 


[ms] 


[/is] 


[/is] 


[moV] 


Y2 


4. 


-0.104 (0) 


4. 


0. 


- 1 




4. 


-0.146 (0) 


4. 


0. 


- 1 


Zr2 


1.984 


-0.062 


2. 


-1.710 (2) 


+ 25 




1.984 


-0.116 


2. 


-1.744 (2) 


+ 76 


Nb2 


1.972 


-0.044 (0) 


2. 


0. 


- 9 




1.972 


-0.080 (0) 


2. 


0. 


- 10 


TC2 


1.924 


0.052 (?) 


2. 


0.012 


- 31 




1.924 


0.084 (?) 


2. 


0.012 


- 34 


Rha 


3.926 


0.630 


3.978 


2.116 ({2}) 


+ 104 




3.928 


1.088 


3.976 


2.122 ({2}) 


+ 196 


Pd2 


1.980 


0.532 (0) 


1.938 


0.862 


- 5 




1.984 


0.762 


1.940 


1.106 (0) 


+ 10 



